function f=R(x) 
f=-(1.078163*x(1)+1.092909*x(2)+1.119773*x(3)+1.123636*x(4)+1.121318*x(5)+1.091773*x(6)+1.141227*x(7)+1.128955*x(8));
